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An algorithm is proposed which transfers the quantum information of a wave 

function (analogue signal) into a register of qubits (digital signal) such that n qubits 

describe the amplitudes and phases of 2 n points of a sufficiently smooth wave func- 

\ tion. We assume that the continuous degree of freedom couples to one or more 

■ qubits of a quantum register via a Jaynes Cummings Hamiltonian and that we have 

O ' 

I/-) . universal quantum computation capabilities on the register as well as the possibility 



o 

a: 



to perform bang-bang control on the qubits. The transfer of information is mainly 
based on the application of the quantum phase-estimation algorithm in both direc- 
tions. Here, the running time increases exponentially with the number of qubits. 

3 : 

Qi. We pose it as an open question which interactions would allow polynomial running 



X 



time. One example would be interactions which enable squeezing operations. 

I. INTRODUCTION 

Traditionally quantum computing and quantum cryptography have been formulated in 
a digital setting, i.e. with qubits [lj]. However, also models with continuous variables have 
been proposed |2| . Protocols for continuous variable cryptography have been investigated 



in detail (e.g. |3]). In ^| several operations on hypothetical continuous variable quantum 
computers have been proposed which can generate arbitrary unitaries. The author argues 
that continuous models possess various advantages compared to the standard model quantum 
computer. Therefore an interface between continuous and discrete registers is desirable, 
since with this device one could combine the advantages of both approaches. There are 
also other reasons why the bridge between continuous and discrete degrees of freedom is 
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an interesting issue of research: The possibility to transfer the wave function of a massive 
particle or the state of a light mode to a quantum register would allow to use algorithmic 
measurement schemes like those proposed in [^] for POVM measurements on the continuous 
degree of freedom. Similarly, the ability to transfer quantum information from digital to 
analogue would allow to use state preparation algorithms in quantum computers p] for 
algorithmic state generation in the analogue system. Furthermore, the implementation of 
POVM measurements with an uncountable number of outcomes on a finite dimensional 
system is only possible if one couples it to a continuous degree of freedom [7 1 . 

An interesting system where the state of a light field is transfered to the state of many 
two- level systems and vice versa is the micromaser (see [s| and references therein). The 
two-level atoms cross a cavity one after another such that at most one atom is present in 
the cavity at any time. While it is passing the cavity, each atom is interacting with the 
cavity field mode via a Jaynes Cummings Hamiltonian. One can prove |8j that every state 
of the field mode can asymptotically be prepared as a limit if an infinite number of atoms, 
initialized to an appropriate state, passes the cavity. It has been shown that for many 
interesting examples small numbers of atoms are sufficient to prepare the desired state with 
high fidelity. Since the final state of the field in the asymptotic scheme does not depend on 
its initial state, the latter has been completely transferred to the outgoing atoms. Therefore 
the system realizes asymptotically the transfer of quantum information in both directions. 
However, the fact that these statements refers to asymptotic behaviour indicates already 
that the state is typically not encoded on a minimal number of atoms. 

Another system where quantum state transfer between a multi-photon state and the 
states of atoms and vice versa has already been experimentally implemented is described 
in j^. In this "quantum memory for light" the eigenvalues of the total spin operator J z 
of the atoms define the basis states of the atomic memory. However, this scheme encodes 
an n-photon states in a collective polarization of an atom ensemble where the number of 
atoms is also of the order n, i.e., the number of qubits is in the order of the dimension of the 
encoded space. In this article we propose a quantum analogue-to-digital converter where the 
number of qubits needed grows only logarithmically in the dimension of the encoded space 
for the cost of an exponential running time of the algorithm. We shall discuss later whether 
this shortcoming can be removed. The question of the cost of accurate A/D conversion is 
directly connected with the question of the computational power of analogue computers, 



which is already an interesting problem in classical computer science [10J. Whether or not a 
continuous degree of freedom could be used to store a "reasonable number" of qubits depends 
on the ability to access a subspace of exponentially large dimension on a "reasonable" time 
scale. 

Here the continuous degree of freedom is represented mathematically by the Hilbert space 
L 2 (IR), the set of square integrable functions on the real line. It is isomorphic to Z 2 (No), the 
space of square summable sequences over No by choosing the eigenf unctions of a harmonic 
oscillator as complete orthogonal system. This shows that continuity or discreteness is here 
not a property of the Hilbert spaces but rather of the considered observables. Our A/D and 
D/A converters refer explicitly to a discretization with respect to a a variable with continuous 
spectrum, e.g. the position variable of a Schrodinger particle, but also applies to the formally 
equivalent variables of a light mode. Using the above isomorphism Z 2 (No) = L 2 (R), it would 
be straightforward to transfer the information such that the state with j oscillation quanta 
is mapped onto the jth binary word in the discrete register. However, here we would like 
to represent the values of the wave functions at 2 n points directly by the coefficients of 
the binary words of the discrete register. For doing so, we restrict us to Schrodinger wave 
functions that are contained in the interval [0, L] (except e. g. exponential tails). The n- 
qubit register is represented by (C 2 ) n with basis states \j) with j = 0, 1, . . . 2 n — 1. Then we 
demand that every sufficiently smooth wave function x i— > ip(x) is converted to the quantum 
register state 

3=0 

in an approximative sense. We will show below that "sufficiently smooth" means that the 
L 2 norm of the derivative of the wave function is not too large. The conversion operations 
that we will use are unitary transformations on 

Now we describe the model in which conversion from analogue to digital is possible. This 
provides us with the available resources for the conversion algorithm. 



We assume that the interaction is described by the Jaynes Cummings Hamiltonian 11 
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= c'^(^W + <7?®a) , (2) 

3=0 



4 



where d > is the interaction strength and we have used the conventions 

a:=±{x + ip), J:=±=(x-ip), a£ := i ( <jf ± z^) ) , (3) 

where <7«' ) denotes the Pauli matrix a a acting on qubit j and x and p are the position and 
momentum operators, respectively, defined by 

(xtp)(x) := xip(x) , (p4>)( x ) '■= —i -r4>( x ) 

ax 

We can rewrite the Hamiltonian in eq. (j2J) as 

n-l 

H = c^2(a ( J ) ®x - a® ® p) . (4) 
i=o 

Note that we choose the oscillator parameters such that mu = 1 and set h = 1 throughout 
the paper. The Hamiltonian (j2J) appears often in physical systems when the continuous 
degree of freedom is an harmonic oscillator, e.g., an oscillation mode of ions in a trap 
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To achieve A/D conversion it is not sufficient to use just Hamiltonian evolution with 
Hamiltonian (J3J), we will also need various other unitary operators. For example, below we 
want to use the terms c4 <8> x and <g> p of eq. (@J) separately. Fortunately, there exists 
already a well-developed technique which allows to simulate various effective Hamiltonians 
QjQl- Propagating the system only for short time intervals with the Hamiltonian (j3J) and 
interrupting this by one qubit unitaries, we can cancel or modify terms of the Hamiltonian. 
We use the fast control limit (also called bang-bang control), i. e. Hamiltonian evolution is 
neglected during one qubit operations are applied. In section 11111 we will explicitly outline 
the one qubit operations and pulse sequences that entail the desired modifications of the 
Hamiltonian. Finally, as our last resource we assume that on the quantum register, we have 
the ability of universal power of quantum computation. Even though we use a specific inter- 
action between continuous and discrete register as a resource for the conversion algorithm 
there are several generalizations which will be obvious after having discussed our method. 
First, the particle wave needs not necessarily interact with all qubits simultaneously and 
with the same strength, one could also have different coefficients. We will furthermore see 
that the only requirement is that one of the interactions a z ®x and a z ®p can be simulated, 
because the other can be obtained by implementing a Fourier transform to the continuous 
system. 
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We now describe the organization of this article. In section |H] we explain the algorithm 
that achieves the conversion of quantum information. Each step of the algorithm is given 
with its corresponding operator that acts on the tensor space of qubit register and wave 
function. In section IHII we describe the procedures for simulation of Hamiltonians which 
generate the required effective Hamiltonians from the given one. In section HVl we describe 
briefly that the time reversed implementation can in principle be used for a digital to ana- 
logue converter. We summarize and discuss our results in section The appendix gives a 
proof of eq. (J3TJ). 



II. THE A/D CONVERSION ALGORITHM 

We first sketch the general idea of the function ing of the A/D converter. It uses a variant 
of the standard phase estimation algorithm Q, Q] in order to bring the wave function 
amplitudes ip(x) into the appropriate place of the qubit register (cf. the scheme of eq. ((H)). 
The essential principle is that the interaction a z <g> x implements a controlled-exp(— ixT) 
operation which allows to use the qubit register als "measurement apparatus" for x. After 
this procedure the joint quantum state displays a high degree of entanglement between its 
qubit and its wave function part. Therefore, in a final step we displace - depending on the 
value of the qubit register - all parts of the wave function to the same location, so that all 
quantum information is deleted in the continuous Hilbert space and transferred to the qubit 
register (again, in an approximate sense). The controlled displacement is done by a a z <S> p 
interaction. 

Before we start the conversion process the phonon wave function ip(x) is contained in 
the interval [-L/2, L/2}. Here the length L should be estimated in such a way that the 
substantial part of the wave function is contained in this interval. We start with the following 
product state 

|^W>:=|1,1,...1>®I$ • (5) 

To make subsequent procedures simpler we displace the wave function sm amount of 
L/2 to the right such that the new wave function lies in the interval [0, L\. The displacement 
operator that achieves this is 

exp ( - % ) . (6) 
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As will be recalled in section ITTT1 we can cancel unwanted terms in eq. (jlj) by standard decou- 
pling techniques by interspersing the natural evolution with one qubit control operations. 

n-1 

D T := exp {iTcJ2 a z } ® P) ( 7 ) 

The application of this operator to the joint state \4>^) (cf. eq. (J5J)) of qubits and wave 
function yields 



) :=D T (|1, $)) = |1, ® exp (-inTcp) \i>) . (8) 



A comparison with formula (JOJ) shows that by choosing the time span T = L/(2 n c) we can 
realize the desired displacement of the wave function. The displaced state is denoted by \ip). 

In the first part of the phase estimation algorithm the n qubits are in a uniform superpo- 
sition of computational basis states and control the application of the operator exp(27u x/L) 
to the wave function. This can be formulated as 

u ( ^ E ^ ® w ) . ( 9 ) 

fc=0 

where operator U is given as 

n-1 

U := exp fern ^ 2 j P 3 ® ^) . (10) 

j=0 

Here Pj is the projection operator that acts on the jth qubit of the qubit register defined 
by 

P j = l®l---®\1)(1\®---®1 . (11) 

Unfortunately, we cannot implement the operator U as it is with our available resources. 
But since P in eq. (fTTjl can be written as P := (1 — er 2 )/2, we can split U into two factors 

n-1 

U = exp 7iiJ2 23 °z ] ®t) exp(7ri(2 n -l)^) =: U R . (12) 

3=0 

The second factor R acts only on the wave function multiplying it with a position-dependent 
phase. Therefore eq. Q is equal to 

2 n -l 

^ ( 2^ E l*> ® * W ) > (13) 

fc=0 

We can realize this transformation in the following three steps 
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i) We first multiply the wave function i/j(x) with the phase 



x 



exp (vrz(2 n - 1) - ) , x e [0, L) , (14) 



which is done as follows. In section ITTT1 we will demonstrate how to realize the operator 

n-1 

D T = exp(-iTc^V<, i) ® x) (15) 

3=0 

with our available resources. Acting with Dt on the joint quantum state |</>^) 
(cf. eq. dHJ)) achieves the phase multiplication 

D T (|1, 1,...1)®|^» = |1, exp (inTcx) |V) • (16) 

ii) To bring the qubit register into the uniform superposition of all computational basis 
states \k) we apply to each qubit the operator 




since 

n-1 _ 2™-l 



E 



h !,...!) = -L \k) • (17) 



2 n l' 

j=0 k=0 

iii) The structure of the operator U as defined in eq. ()12|) is very similar to the one of the 
operator Dt of eq. 1)1 5|) . The only difference are the factors 2 J in the exponent of U. 
Clearly we cannot increase the strength of the interaction by any selective decoupling 
scheme. In order to obtain a unitary which would correspond to the exponentially 
growing interaction we need exponential interaction time (see section |TTTJ). 



After the application of these steps the quantum state of the joint system is changed to 

on 1 r 

V- f , f,,. 1 /27rikx^ 



■ i2) ) = f/|0 (1) ) = E / ^{l fc >®2^ ex P (^r^) ^) 1^) } • ( 18 ) 

fc=0 ^° 



Note that we have preferred to use the notation \x) even though position eigenstates do not 
exist (readers who appreciate mathematical rigor may forgive us). The whole expression is 
nevertheless a well-defined state in the joint Hilbert space. The second part of the phase 
estimation algorithm consists of the application of an inverse Fourier transform JF _1 to the 
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qubit register. The Fourier transformation (and its inverse) can be efficiently implemented 
on a quantum computer Q, Q] . After this transformation the quantum state becomes 



10(3)) = (r 1 ®!) |0( 2 )) 



2 n -l 



1=0 



2"-l 



E / ^{^)io®^E 



exp 



fc=0 



27T2 k 



\x) 



(19) 



|</>^) displays a high degree of entanglement between its qubit and its wave function part. 
The following procedure removes a large part of this entanglement and thus completes the 
transfer of quantum information. For this purpose we apply the following operator to the 
quantum state 



V = exp 



iL 



n-l 



(20) 



3=0 



As before we will discuss the implementation of this operator with our resources in section 
TTT1 Acting with operator V on the quantum state of eq. (fT9*j) . the wave function part 
is displaced where the amount depends on the entangled qubit state |/) 

2" — 1 l 

) : = V|0 (3) > := E / dx\i>{x) |0 <g> 

,_n JO 



«=0 



2™-l 



exp 



fc=0 



where we have used the quantity 



2 77/ /,:[ - - — 

L 2™ 



h := L 



\x + h-{l L)/2 n ) 



(21) 



2 n+l 



(22) 



Making the substitution x' := x — (I L)/2 n we can rewrite this quantum state as 



1=0 



where 



2 " 1 1 1 1- exp (27rc2"y) 

^) : = E 2^ exp(2*,*y) = 1 _^ {2my) 

The function is periodic with period 1. One can easily show that 



dy\g{y)f 



(23) 



(24) 



(25) 



Besides, the function \g{y)\ 2 becomes highly peaked around y = for a large number of 
qubits as shown in fig. [T] In fact, the width of \g(y)\ 2 is proportional to 2 _n . Outside the 
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FIG. 1: The function |<?(:r)| 2 for n = 6. 

peak the function g(y) takes on values that have only a small modulus (oc 2~ n / 2 ). Thus the 
wave function part in the quantum state \4> <y4}> ) of eq. (|23|) displays a peak around x' = h 
where h (cf. its definition in eq. ([22))) is approximately the mid point of the interval [0, L]. 

The result in eq. (|23|) is almost satisfactory, but a minor technical point should be men- 
tioned. We would like that in eq. (|2~Hj) each qubit register state \l) has the same "wave 
function" in its corresponding continuous part of the tensor space. We have seen that the 
function g(x' / L) can be neglected outside its peaks. Considering the /-dependent integra- 
tion bounds in eq. (|23|). the relevant wave function is not the full peak of g[x' j L) only for low 
and high values of I E {0, 1, . . . , 2 n — 1} . Hence we have to make the additional assumption 
for the original wave function that ip(x) ~ holds in the two subintervals of length O (2~ n ) 
that join the end points and 1 of the interval [0, L\. 

Using this technical assumption and the fact that g(x'/L) can be neglected far away from 
its peak, we can rewrite the quantum state of eq. (|23|) to a very good approximation as 

2 n 1 /•W' i 7 r / 

m*J2 dx '{^(^+^\l) ®9{j)\h + J)}, (26) 

1=0 J ~ w 

where W > is some multiple of the width of the function g(x'/L) and hence W oc 2~ n . 
Eq. (|2fj)l shows that now the quantum information has been transferred to the qubit register, 
since the continuous Hilbert space is left with the "standard wave function" g(x'/L). This 
is in accordance with the no cloning theorem that precludes the copying of quantum 
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information. 

The result eq. ()26|) still shows some degree of entanglement as qubit and wave function 
part are connected via the integration over x' . In order to assess the magnitude of this 
entanglement, we calculate the reduced density operator of the qubit register 

/oo 
dx (x|0 (4) ) (0 (4) |x) 
-oo 

= f E />{KDN^H^+*)} • < 27 > 

In the following we show that for a large number of qubits this density operator can be 
replaced by the density operator 



n dr sr-^ (i L\ - (k L 

P ... 

which is the density operator of a pure state |\&) with 



The constant d = 1 + 0(2~ n ) ensures the correct normalization in eqs. (|28J) and (|29|) . To 
demonstrate the possibility of replacing p by p°, it is appropriate to consider the trace norm 
I22I ] of the difference operator Ap := p — p°. For - thanks to Holder's inequality - we can 
bound expectation values for an observable A as 

Tr[(p-p°)i]| < \\p-p% ||i|U , (30) 

where || ||oo is the operator or spectral norm of A. In the appendix to this article we show 
the following bound for the trace norm 

\\p-p°\\i<b^, n»l, (31) 

where 

/1 1 
dx \tp'(x)\ 2 y . (32) 

Eqs. (I30|) and (J31|) show that in the limit of a large number of qubits the density operators 
p and p° are equivalent for all observables A whose || H2 norm diverges slower than 2 n . 
We thus need a large number of qubits n in order to represent the quantum information 
of a wave function faithfully in a qubit register. Furthermore, according to eq. (}3*2*j) the 
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accuracy of this representaion is also determined by the L 2 norm of the derivative of the 
wave function. Hence the smaller the derivative of the wave function, the better works its 
conversion into digital information (for a fixed number of qubits). 

We have already mentioned that the time needed for the execution of the phase estimation 
algorithm grows exponentially with the number of qubits. Thus there is a trade-off between 
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e. op- 



accuracy and speed for our A/D conversion algorithm. Squeezing operations 
erators of the form 

S(r):=exp(^[a 2 -(at) 2 ]), rER, (33) 

could speed up the phase estimation algorithm considerably. For, the squeezing operator 
could magnify the wave function %p(x) by a factor A > 1 while preserving its shape. This 
would decrease the time for phase estimation by a factor A -1 . Note that we can rewrite S(r) 
as 

S(r) = exp ir(xp + px) . (34) 
In order to generate such unitaries we have to simulate the Hamiltonian 

xp + px . 

To achieve this, we observe that 

[<j x <8 x, a y <g> p] = ^a z <g> (xp + px) . 

We conclude that (xp + px)/2 can be obtained by the following second-order simulation 
scheme which applies the following 4 Hamiltonians for a small time AT: 

(l)a x ®x, (2)a y ®p, (3) -a x ®x, (4) -a y ®p . (35) 

Up to terms 0((AT) 3 ) we obtain a time evolution according to the desired Hamiltonian 
multiplied with a slow-down factor (AT) 2 provided that the qubit is set to the state |1). 
Due to 

(1 ® S(r) k )(a z <g> x)(l <g> S(~r) k ) = a z ®r k x 

one could simulate exponentially large interaction time by a linear number of concatenated 
squeezing operations before the interaction has taken place and undoing the squeezing af- 
terwords. However, the problem with a second-order simulation is that the running time 
increases with the desired accuracy. Since the required error decreases exponentially with 
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the desired qubits we expect here also exponential running time. However, on a scale where 
squeezing operations are available with sufficient accuracy one could nevertheless expect a 
speed up. 



III. SELECTIVE DECOUPLING AND SIMULATION OF HAMILTON! ANS 

Simulation of Hamiltonians by interspersin g th e natural time evolution with fast con- 
trol operations is used in NMR since decades [l_2| . These techniques are subject of many 
theoretical investigations 2^3,0]. Here we refer only to very basic ideas. 

Let H be the natural Hamiltonian (J2J). Using the anti commutator relation between Pauli 
matrices 

{di, crj} = , j e {x, y, z] , (36) 

we get the equation 

n— 1 n— 1 

(gK exp ( - i AT H) (gK = exp ( - i AT H') , (37) 

j=0 j=0 

where 

n-1 

H':=cJ2{- °x ] ® £ - ^ ® P) • (38) 
j=o 

The operators H and if' do not commute, but for a small time interval AT < 1 we can use 
the Baker Campbell Hausdorff formula 

exp ( - i AT H) exp ( - i AT H') = exp ( — i AT (H + H') + 0(AT 2 ) ) , (39) 

where 

n-1 

H + H' = -2 c a y ] ® P ■ ( 4 °) 

j=0 

Thus to leading order in AT the unwanted ax <8> ^ terms have canceled each other. 

Therefore, if during a time interval T we change frequently between Hamiltonian evolution 
and the product of one qubit operations ®™zla y , we can realize the unitary operator 

n-1 

B T = exp(iTcJ2 a y ] ® P) • ( 41 ) 

j=0 
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In the language of |21| we have now "simulated the Hamiltonian" 

c yj a® ® p . 

j=0 

In a similar way, we can also select the term with x in (j2J by applying cr x .-rotations to all 
qubits. Complete decoupling can be achieved if we apply a z to all spins since this reverses 
the sign of the x and the p term. If we want to cancel all terms except from the interaction 

a® ® x (42) 

for one specific qubit j we apply a x to qubit j and a z to all the other qubits. To simulate 
the time evolution 



exp 



3 

we may concatenate the commuting unitaries 



exp (iT2^a^' (g> x 



which are obtained by applying the simulated interaction (|42j) on a time interval of length 
2 J 'T. If we apply arbitrary single qubit unitaries to qubit j initially and apply Uj after- 
words, we obtain the time evolution 

exp (iT2 j {U j a < j ) U}) ® £ 

This shows that we can replace the Pauli matrix a x in eq. (|42|) by other Pauli matrices or 
by — <y x as we like. 



IV. DIGITAL-ANALOGUE CONVERSION 



Using the results of the previous sections, we can easily describe an algorithm for digital- 
analogue conversion. Roughly speaking, the argument is as follows. Let Z be the transfor- 
mation on the continuous and the discrete degrees of freedom which implements the complete 
analogue-digital conversion algorithm. Provided that the wave function ip was sufficiently 
smooth the system ends up almost in the product state 



(43) 
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where is a superposition state with coefficients L/2 n ) with j = 0, 1, . . . ,2 n — 1 
as defined in eq. (J29|) and ip is the function whose absolute square was plotted in fig. [T] 
translated by L/2. If we apply to the state (|4*H|) we obtain hence almost the original wave 
function ip. It is clearly required that the values L/2 n ) correspond to some sufficiently 
smooth wave function. We have explained that sufficiently smooth means here that |^'(x)| 
is small enough. Hence the vector corresponds to a smooth wave function whenever the 
values L/2 n ) do not vary too much for adjacent j. 

In order to obtain definite accuracy bounds from this idea we recall that the reduced 
state p of the discrete register satisfies 

l|p-|*X*ll|i<6^=:c- (44) 
If we set A = \SS?)(SSf\ in eq. (|30j), we obtain 

||p-|*)(*|||i>|W|*)-l|. (45) 
Hence the largest eigenvalue of p is at least 1 — e. Let 

\4>) = Y,Ci\<*i)®\Pi) (46) 
j 

be the Schmidt decomposition of the exact bipartite state after the A/D conversion. Then 
the absolute square |c | 2 of the dominating coefficient c is the largest eigenvalue of the 
reduced density operators on both subsystems. Hence the square of the norm distance 
between \(j>) and co|«o) <8> |/3o) is at most e. It follows that the error which arises from 
replacing the joint state by the tensor product state 

|ao)<ao|® IA))(/3o| (47) 

is of the order e. We know furthermore that we may replace the state of the digital system 
by such that the error in trace norm is in the order of e. We conclude that applying W 
to 

|a )(a |®|*)(*| (48) 

leads to the wave function ip up to a trace norm error in the order of e. The initialization of 
the continuous system to the state \ocq) is clearly a non-trivial task. It is some wave packet 
which is similar to the function g translated by L/2. To construct algorithms which work 
also for more general initializations shall not be our subject here. 
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V. CONCLUSIONS 

In this article a quantum algorithm has been outlined that can read in the quantum infor- 
mation of a wave function into a qubit register. This can be viewed as the quantum analogue 
of an A/D converter which encodes an A-dimensional system obtained by discretization of 
a continuous wave function into a log 2 N qubit register. 

The principal ingredient for the A/D converter is the use of the phase estimation al- 
gorithm that is already popular for other purposes in quantum computing. The principal 
resources are interaction Hamiltonians which are tensor products of a position or a momen- 
tum operator with a Pauli matrix of the discrete system. Effective Hamiltonians of this 
type can for instance be obtained from the Jaynes Cummings Hamiltonian using standard 
techniques for selective decoupling. 

Whether the proposed A/D and D/A converters could be realized depends on the one 
hand on the experimental progress. On the other hand and perhaps even more important, 
further theoretical work is required to investigate more systematically the advantages that 
arise when one combines discrete and continuous degrees of freedom in quantum computing. 
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APPENDIX 

In this appendix we prove eq. (|3ip. To simplify the notation we set L — 1. We start with 
the density operator p and use the mean value theorem in eq. P7|) 

^ a + x ) = ^ (io + x ^ a + ^ x) ) with ° - x) 1 - ixi - w (49) 
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elements of the difference operator Ap = p — p° in the standard basis as 

l /"' . r ,, r / / \ - , / k 



and analogously for the conjugate function ip(k/2 n + x). Therewith we can_write the matrix 

-w 
-w 



( A P) 3 k = ¥ / dz{\g(x)\*xty(±) + xj) 



(50) 



We can rewrite this expression in an operator form 

1 f w r r n 

Ap = ¥ J w dx {\9(x)\ 2 xp)(K\ + \K)(*\\} , (si) 

where the 2 n dimensional vectors | and | ty' x ) are defined in the standard basis as 
Using the triangle inequality, eq. (|5Tj) leads to the following estimate for the trace norm 



i / r w 

||Ap||i<-( / dx\g 
z v J-w 



x)\ 2 \x\ 



sup {|| I*) (%\ + \%) (*| ||i} . (53) 



-W<x<W 

We now consider the trace norm of the operator 

J3* := | ¥> | + | tf.) <¥ | • 

We see from this formula that the operator acts non-trivially only on the two dimensional 
subspace span{ | \ ty' x ) }. Choosing an ONB {| ), 1 1 )} for this subspace with | ) := 
I ^) I II I ^) II 2 ( II II 2 is the Euclidean norm), we arrive at the following matrix representation 
for B x 

f (e^ | ) + < | e x ) (e*|l) 
(lie,) 

where | e x ) := | ^)/|| | ^' x ) |b is a unit vector. From eq. (j54j) the trace norm of B x can be 
bounded uniformly 

l 

ii b x Ik < ii i *> ib ii i *;> || 2 \ B M ^ 4 ii i *> n= ii i ^) h ■ (55) 

j, fc=0 

In addition, for n > 1 we can evaluate the Euclidean norms of the vectors | \&) and | ^f' x ) 
approximately as 



F e = |||*>|b |||^> lb 



(54) 



*) 



2~ 



22 



dx \ip'(x)\' 



(56) 
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Using the results of eqs. (|B3j) and (|5l)|). we can bound the difference between the two density 
operators as 

^dx \g(x)\ 2 \x\) ( J dxW{x)\ 2 y . (57) 

Since for n 3> 1 the width W oc 2 _ri and thus 

-w 



dx {\g{x)\ 2 \x\} oc2~ n , (58) 



li- 



the inequality (jBTjl establishes the bound eq. (j^Tj) with the constant 



6oc ( / dx \^'{x)\ 2 f . (59) 
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